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Abstract 

We study cosmological implications of supersymmetric axion models in which the ax- 
ion scale is generated radiatively. Such models lead to the so-called thermal inflation and 
subsequent reheating should be constrained not to yield a too large axion energy density 
at the time of nucleosynthesis. We examine how plausible it is that this nucleosynthe- 
sis constraint is satisfied for both hadronic and Dine-Fischler-Srednicki-Zhitnitskii type 
axion models. Baryogenesis and the possibility for raising up the cosmological upper 
bound on the axion scale in thermal inflation scenario are also discussed. 



One of the attractive solutions to the strong CP problem is to introduce an anomalous 
Peccei-Quinn (PQ) symmetry U(1)pq |IJ. This solution predicts a pseudo-Goldstone boson, 
the (invisible) axion R [| , whose decay constant F a is tightly constrained by astrophysical and 
cosmological arguments. The allowed band of the axion scale F a lies between 10 10 GeV and 
10 12 GeV |4[] which is far away from the already known two mass scales, the electroweak scale 
and the Planck scale Mp = 1/ '\/WkG^. It is certainly desirable that this intermediate scale 
appears as a dynamical consequence when the known mass scales are set up in the theory. 
This indeed happens || in some class of spontaneously broken supergravity models which 
are commonly considered as the underlying structure of the supersymmetric standard model. 
Such models typically contain two basic mass scales, Mp and the scale of local supersymmetry 
breaking Ms in the hidden sector leading to 7773/2 = Mj/Mp = 10 2 ~ 10 3 GeV. Supergravity 
interactions then generate soft supersymmetry breaking terms in the supersymmetric standard 
model sector which are of order 777,3/2. In this scenario, radiative corrections to the Higgs 
doublet mass-squared associated with the large top quark Yukawa coupling can naturally lead 
to the electroweak symmetry breaking at the scale M\y ~ 7773/2- When the PQ fields which 
are responsible for the spontaneous violation of U(1)pq correspond to flat directions of the 
model, the intermediate axion scale F a can also be radiatively generated in terms of Mp and 
7773/2. In such a scheme, as was recently emphasized, the early universe experiences the so- 
called thermal inflation and subsequently a period dominated by coherently oscillating flaton 
fields [[J. The aim of this paper is to examine cosmological implications of PQ flatons in 
supergravity models with a radiative mechanism generating the axion scale. 

One possible cosmological consequence of PQ flatons is the impact on the big-bang nucle- 
osynthesis through their decay into axions. In the scheme under consideration, PQ flatons have 
generally order-one coupling to the Goldstone boson (the axion) in the unit of 1/F a 0. As we 
will argue later, axions produced by decaying flatons are hardly thermalized. In this paper, we 
first consider the energy density of these unthermalized axions at the time of nucleosynthesis 
together with its implications for both hadronic axion models and Dine-Fischler-Srednicki- 
Zhitnitskii (DFSZ) type models 0. Even when one takes a rather conservative limit on the 
axion energy density, this consideration provides a meaningful restriction for generic hadronic 
axion models and also for DFSZ type models with a rather large flaton mass. 

As another cosmological implications of PQ flatons, we consider the possibility of rais- 
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ing up the cosmological upper bound on the axion scale F a through the late time entropy 
production || by oscillating PQ flatons. We argue that F a can be pushed up to about 10 15 
GeV without any cosmological difficulty in thermal inflation scenario. Finally, we point out 
that the Dimopoulos-Hall (DH) mechanism |J for a late time baryogenesis can be naturally 
implemented in thermal inflation scenario. In the conclusion, we note that the case of n = 2 
or 3 [see Eq. (3) below] provides a very concordant cosmological scenario. 

We begin by describing how the intermediate axion scale can be radiatively generated in 
supergravity models in which the PQ fields correspond to flat directions. Let us consider a 
variant of the model of Ref. with superpotential 

W = k^l + + h N NN(pi + h L LH 2 N + -.. (1) 

where H\p, are the usual Higgs doublets, N is the right-handed neutrino component and the 
ellipsis denotes the supersymmetric standard model part of the superpotential. In order to 
implement the PQ symmetry, two gauge singlet superfields <pi 2 with PQ charges qi t 2 are 
introduced. The structure of the superpotential is determined by the PQ charge assignment: 
q 2 = —(n+2)qi, q H +q H = —(n+l)qi and so on. Obviously the PQ fields 4>\ and 2 correspond 
to flat directions when nonrenormalizable interactions and supersymmetry breaking effects are 
ignored. This model can be considered as a supersymmetric generalization of the DFSZ axion 
model (but endowed with a radiative mechanism generating the axion scale) in the sense that 
the Higgs doublets carry nonzero PQ charges. Note that the second term in the superpotential 
yields the correct scale for the Higgs mass parameter \x = h{<pi) n+l / Mp upon spontaneous 
breaking of the PQ symmetry [K|. Taking into account the radiative effects of the strong 
Yukawa coupling h^NNcpi, the soft mass-squared of <pi becomes negative at scales around 
F a ~ (0i ), and thereby driving 0i to develop vacuum expectation value at an intermediate 
scale. This Yukawa coupling is also necessary to keep the field (pi in thermal equilibrium at 
high temperature T > mi for which (0i) = 0. Neglecting the field 02, the renormalization 

group improved scalar potential for the singlets is given by 

Iflii 1 2,1+4 

V = V -ml\^ + k^-^_^ (2) 

where m\ is positive and of order rriy 2 , and Vo is a constant of order m^^F^ which is introduced 
to make V((<j)i)) = 0. Clearly the minimum of this scalar potential breaks U(1)pq by 

(0i> ~ F a ~ (m 3/2 M^ n +\ (3) 



where we have ignored the coefficients of order unity. The integer n fixes the size of the axion 
scale. For the smallest value n — 1, the axion scale F a ~ yJmy^Mp fits into the usual allowed 
band of the axion scale: 10 10 GeV ~ F a ~ 10 12 GeV. Later we will argue that the upper bound 
on F a can be relaxed and thus a bigger value of n is allowed also. 

The above radiative mechanism generating the axion scale has substantial influence on 
the history of the universe || [□]]. At high temperature, 0i receives a thermal mass b~m\ ~ 
|/itv| 2 T 2 ^> m\ leading to (0i) = 0. This thermal mass is generated by right-handed neutrinos 
in the thermal bath. Note that the right-handed neutrino N becomes massless when (0i) = 
and thus copiously produced when T ^> mi. During this period, (02) = also. When the 
temperature falls below T ~ V^ 4 , which is about ^Tn^/ 2 F a , the universe is dominated by 
the vacuum energy density Vq and thus there appears a short period of thermal inflation. 
Below T < mi ~ m 3 / 2 , the effective mass of 0i becomes negative and then X develops an 
intermediate scale VEV given by Eq. (0). With (0i) — F a , the other flaton field 02 develops 
also a VEV of order F a through the A-type soft SUSY breaking term, fcv40™ +2 2 /Mp, in the 
scalar potential. This procedure makes the thermal inflation end and subsequently the early 
universe experiences a period dominated by coherently oscillating PQ flaton fields 0i and 02. 
More precisely, the oscillating flaton corresponds to a combination of the two complex scalar 
fields 0i and 02 which is orthogonal to the axion field a = J2i Ci arg(0j) where q = qi(4>i) 2 /F a . 



NS bound. After the period of coherent oscillation, the universe would be reheated by 
the decay products of the oscillating flaton ip. A feature peculiar to the PQ flatons is that 
their decay products include axion as one of the main components 0, |TJ]. The energy density 
of these axions at the time of nucleosynthesis (NS), (p a )NSi should satisfy the conventional 
nucleosynthesis bound on the extra energy density: 

(^) < *N U . (4) 

V Pv ) NS 

Here p u denotes the energy density of a single species of relativistic neutrino and 5N U is the 
number of extra neutrino species allowed by nucleosynthesis. In the past, 5N U has been argued 



to be 0.3 or even smaller as 0.04 [fLlq| . However, although claimed to be quite conservative, 
more careful recent analyses do not exclude even SN V = 1.5 |Tj|. Here we do not take any 
specific value of 5N U , but examine the implications of the above NS bound for 5N V = 0.1 ~ 1.5 
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with the hope that one can push 5N U down to the value 0.1 in the future. 

Before evaluating (p a ) NS , let us first determine the reheat temperature T RH by parameter- 
izing the width of the flaton decay into thermalizable particles as = B~ 1 M 3 /64ttF 2 . Here 
M 9 denotes the flaton mass and the prefactor B~ l will be presumed to be of order 10, which 
is a proper choice for (p a )NS to satisfy the above NS bound. The reheat temperature is then 
given by 

1 C-l) 1 ' 2 (^L) (^-f GeV, (5) 
yRH V F * \Bj \ F n ) V300GeV/ ' V ; 



RH 



where g RH = g*(T RH ) counts the effective number of relativistic degree of freedom at Trh- 
The entropy production factor S , a fter/>Sbefore for this reheating is of order V /ml^ 2 T RH which 
is roughly of order 10 2 (Mp / 'm 3 / 2 )( 5n ~ 1 ^( 2n+2 \ This huge entropy dumping at relatively late 
time was considered as a promising source for erasing out various unwanted cosmological relics, 



especially, cosmologically dangerous string moduli [I I 



In order to evaluate (p a )vs ; one needs to know whether axions produced by the late flaton 
decay have ever been in thermal equilibrium with the thermalized plasma of normal light 
particles. If axions were in thermal equilibrium at some moment but later frozen out at 
temperature T/, we would have 



'p\ = 4 / 43/4 \ 
,pJns 7 \9*( T f)J 



4/3 



(6) 



However if axions have never been in equilibrium, (p a )NS is simply determined by the effective 
branching ratio B a measuring how large fraction of flatons are converted into axions during 
the reheating. Roughly B a ~ T a /T tot with the decay width r a of (p — > 2a, however axions can 
be produced also by the secondary decays of the decay products of flatons. For unthermalized 
axions at Trh, the ratio between p a and the energy density p r of thermalized radiation would 
be simply B a /{1 — B a ). We then have 

1/3 

(7) 



PA _ 43 B a 4 Z 43/4 
fv) ns 4 1 -- B a7 \g RH J 



In order to see whether axions have ever been in thermal equilibrium, let us consider the 
axion interaction rate T^t = (av)N r where a denotes the cross section for the axion scattering 
off the thermalized radiation with energy density p r and number density N r . A careful look at 
of the reheating process indicates that p r ~ R~ 3 ^ 2 , N r ~ i?" 1 / 2 , and p v ~ R~ 3 e~ tTtot during the 
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reheating period between to and t D ~ where R denotes the scale factor and to corresponds 
to the time when the relativistic particles produced by the flaton decay become the major 
part of the radiation M. A simple dimensional analysis implies that the axion cross section 
can be written as a = ( / ~fi + 7 2 (m/-E l ) 2 )/47r.F 2 , where E denotes the center of mass energy, 7 12 
are dimensionless constants of order unity or less, and m corresponds to the mass of target 
particle. We then have (av) = (7i + 7 2 (m/(i^ )) 2 (i?/i?o) 2 )/47r.F 2 . With the informations given 
above, it is straightforward to see that the ratio Z = T int /H is an increasing function of R 
during the reheating period of t < t < t D . 

Let us now consider the behavior of Z = T int /H after the reheating. At t > t D with 
T < Trh-i the entropy production almost ends and thus N r ~ R~ 3 ~ T 3 , H ~ R~ 2 ~ T 2 as in 
the standard radiation dominated universe with an adiabatic expansion. Using Eq. (3) with 
M v ~ m 3 / 2 and Eq. (5), we find 

^ = 3 x i -%y 2 ™^ 2 



H F a 

Jm,J V 10 2 / V M 



/ T \ / n \ fM \ 3 ( n - 1 )/ 2 ( n + 1 ) 



RH 

\2\ 



where J — + 72 (m/ (E)) x 

For n > 1, the above result for t > t D together with the fact that Z = T int /H is an 
increasing function of R during to < t < t D readily implies that Z 1 and thus axions have 
never been in equilibrium. For the case of n = 1, we need a bit more discussion about the size 
of 7. Obviously at tree level, any nontrivial axion couplings to SU(2) x £7(1) non-singlet fields 
arise as a consequence of SU(2) x U(l) breaking. In other words, tree level axion couplings 
to normal fields are induced by the mixing with the Higgs doublets. As a result, tree level 
axion couplings can be described effectively by dimensionless coupling constants which are of 
order m/F a where m corresponds to the mass of the particle that couples to the axion. This 
means that the energy dependent part of the axion cross section, i.e. the 72-part, is due to 
tree level axion couplings, while the energy independent 71-part is due to the loop-induced 
axion couplings like -^fraG^G^ . As a result, 7 is suppressed either by the loop factor (^2) 2 
or by the relativistic factor (m/ (E)) 2 . Then we can safely take 7 ~ 1, implying Z 1 for the 
case of n = 1 also. 

In the above, we have argued that Z = T int /H 1 and thus axions have never been in 
thermal equilibrium. Then the axion energy density at nucleosynthesis is given by Eq. (7) and 
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the NS bound (4) leads to 



B, 



< 0.24 ( 5 M 



1/3 



(9) 



l-B a ~ \ 1.5 J V 4 3/ 4 / 

The above nucleosynthesis limit on B a depends mildly upon the reheat temperature Trh 
through the factor {g RH / g NS ) 1 ^ where g NS = g*(T NS ) = 43/4, while it is rather sensitive to 



the discordant number 8N V which is presumed here to be in the range 0.1 ~ 1.5 (JT2J, [13 
For Trjj above 0.2 GeV but below the superparticle mass, we have g RH /g NS = 6 ~ 10, while 
g RH /g NS = 1 ~ 3 for Trh < 0.2 GeV. We thus have just a factor two variation of the limit 
when Trh varies from the lowest allowed value 6 MeV |1| to the superparticle mass of order 
100 GeV. In summary, the NS limit (9) indicates that we need to tune the effective branching 
ratio B a to be less than 1/3 ~ 0.02 for 5N V = 0.1 ~ 1.5. 

Implication on flaton couplings. We now discuss the implications of the NS bound (9) 
for generic supersymmetric axion models with a radiative mechanism generating the axion 
scale. Since the models under consideration involve too many unknown free parameters, we 
just examine how plausible it is that the NS limit (9) is satisfied for the unknown parameters 
simply taking their natural values. To proceed, let us write the flaton couplings responsible 
for the flaton decay as 

C v = Cpq + Cssm, (10) 

where Cpq describes the couplings to the fields in PQ sector, while Cssm describes the cou- 
plings to the fields in supersymmetric standard model (SSM) sector. Schematically Cpq is 
given by 



£pq = (My + M 2 y + (M^0 + h.c)) (11) 



where a, if', and (p denote the axion, other flaton, and flatino respectively. In the above, we 
have ignored the model-dependent dimensionless coefficients of each terms which are of order 
unity in general. 

The flaton couplings to SSM fields are more model dependent. In DFSZ type models, 
flaton couplings to the SSM sector are essentially due to the mixing with the Higgs doublets. 
Then flaton couplings can be read off by making the replacement 

H t -+ Vi + (12) 
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where Vi = (Hi), Xi/s are model-dependent coefficients which are generically of order unity. 
Then again schematically 

Cssm = {(M lX x' + M 2X X + h.c.) + M\A^ + Mjzz' + M 2 5 \z\ 2 } , (13) 

where z and z' denote spin zero fields in the SSM, e.g. squarks, sleptons and Higgs, with 
their fermionic partners x an d x' i while (A^ 1 , A) stands for the gauge multiplets which become 
massive due to the Higgs doublets VEVs, i.e W and Z. The order of magnitude estimate of 
the dimensionful coefficients leads to: M x ~ M x , M 2 ~ M 3 ~ M w , Ml ~ M X (A + //cot/5), 
and Mf ~ M 2 + cos 2(3, where M x and M w denote the masses of x and W respectively, 
tan/3 = V2/V1, and again we have ignored the coefficients of order unity. 

As is well known, besides DFSZ type models, there are another interesting class of axion 
models named as hadronic axion models. In hadronic axion models, all SSM fields carry 
vanishing PQ charge and as a result flaton couplings to SSM fields appear as loop effects. As 
an example of supersymmetric hadronic axion model with a mechanism generating the axion 
scale radiatively, let us consider a model with 

m+2 j 

W = + hQQQTfa + ■ ■ ■ , (14) 

where <f>i t 2 are gauge singlet flatons, and Q and Q c stand for additional heavy quark and anti- 
quark superfields. Again the soft mass-squared of <pi becomes negative at scales around F a by 
the radiative corrections involving the strong Yukawa coupling hqQQ c <j)i, thereby generating 
the axion scale as (<f>i) ~ (02) — F a . A peculiar feature of this type of hadronic axion models 
is that at tree level flatons do not couple to SSM fields, while there are nonzero couplings to 
PQ fields as Eq. (11). Flaton couplings to SSM fields are then induced by the loops of Q and 
Q c , yielding 

£ssm = (\G; u G a ^ + a a 7 • 9A a )) , (15) 

where {G a , A a ) denotes the gluon supermultiplet (and possibly other gauge multiplets) and 
again we ignored dimensionless coefficients of order unity. 

It is now easy to notice that, due to the loop suppression in Cssm, most of oscillating 
flatons in hadronic axion models decay first into either axion pairs, or lighter flaton pairs, or 
flatino pairs, as long as the decays are kinematically allowed. Lighter flatons would experience 
similar decay modes, while flatinos decay into axion plus a lighter flatino. Then in the first 
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round of reheating, flatons are converted into either axions or the lightest flatinos. The lightest 
flatinos will eventually decay into SSM particles. Because of kinematical reasons, e.g. the mass 
relation M v > 2M^ and the phase space suppression factor (1 — 4M|/M ( ^) 1 / 2 in the decay 
tp — > 2tp, more than half of the original flatons would be converted into axions, i.e. the effective 
branching ratio B a ~ 1/2, unless the flaton coupling to the lightest flatino is unusually large. 
This is in conflict with the NS limit (9) even for the most conservative choice SN V = 1.5, 
implying that hadronic axion models with a radiative mechanism can be compatible with the 
big-bang nucleosynthesis only when the models are tuned to have an unusually large flaton 
coupling to the lightest flatino. 

In DFSZ type models, flatons have tree level couplings to SSM fields which are of order 
M SSM /F a or M g 2 gM /F a where M SSM collectively denotes the mass parameters in the SSM, e.g. 
M t) M W) fx, A, and so on [see Eq. (13) and the discussions below it]. Thus if M v 3> M SSM , 
the reheating procedure would be similar as that of hadronic axion models and then the NS 
limits provides a meaningful constraint on the flaton couplings to the PQ sector. For the case 
that M v is comparable to M ggM , the most conservative choice of SN V = 1.5 would not provide 
any meaningful restriction on DFSZ type models. However it is still nontrivial to achieve B a 
significantly smaller than 1/10. A careful examination of the flaton couplings in DFSZ type 
models suggests that, among the decays into SSM particles, the decay channels to the top (t) 
and/or stop (t) pairs are most important. Flaton coupling to the top (stop) is of order M t /F a 
(M~/F a ), while the coupling to the axion is of order M^/F a . As a result, B a significantly 
smaller than 1/10 implies that the flaton couplings to the top and/or stop are unusually large 
in view of the relation > 2M t (2Mf). One of the efficient way to achieve such a small 
B a is to assume that there is a sort of mass hierarchy between the lighter stop mass-squared 
M~ and the heavier stop-mass squared M~^, allowing for instance 4M~ < M 2 < \M~ . This 
would be the case when Mf + M? ~ M~ c + ~ M t (A + \± cot j3) where Mf and Mf c denote 
the soft squark masses. Since the flaton couplings to stops are determined not only by the 
mass parameters (e.g. M t and A) but also by additional dimensionless parameters x^- defined 
in Eq. (12), the flaton coupling to the lighter stop t\ would be of order M~ 2 /F a , not the 
order of M^jF a . To be more explicit, let us write this coupling as xiM~^ip\ti\ 2 / F a . With the 
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flaton-axion coupling given by x a M 2 {pa 2 /2F a , we find 

r„ 1 f Xa \U2MX(, im?\- 1/2 



< , , , 1 £ . (16) 

h 32 \xj \M h ) V Ml) K J 

This shows that B a can be smaller than about 10~ 2 for the parameter range: x a ~ X\ and 
AM? < M 2 < \Ml . 

tl V 4 t2 

Relaxation of the bound on F a . The reheat temperature can not be arbitrarily low in 
order to be compatible with the big bang nucleosynthesis. Since flatons produce large number 



of hardrons, the bound Trh > 6 MeV has to be B met JL4]. With Eq. (5), this leads to the 
upper bound: 

K , 2xl0 »(°A) l/2 (M^f GeV . (17) 

\B a J V300GeV/ V ; 

Once one uses the relation F a ~ (M !fi Mp) 1 ^ n+1 , this means that only n — 1, 2, and 3 are 
allowed by the big-bang nucleosynthesis. 

As is well known, another upper bound on the axion scale can be derived by requiring 
that the coherent axion energy density produced by an initial misalignment should not exceed 
the critical density [0. If there is no entropy production after the axion start to oscillate 



at around T ~ 1 GeV, this lead to the usual bound: F a ~ 10 12 GeV. When n = 2 or 3, the 
corresponding axion scale F a ~ (M v ,Mp) 1 / n+1 would exceed this bound. However in this case, 
the reheat temperature @ goes below 1 GeV. Then the coherent axions may be significantly 
diluted by the entropy dumped from flaton decays, thereby allowing F a much bigger than 10 12 
GeV!. 

Axion production in matter-dominated universe, e.g. flaton oscillation dominated universe, 
has been considered in Ref. |jl4| , |T6f assuming m a (T) oc T -4 . For our computation, we take 



the power-law fit of the temperature dependent axion mass [ITT 



m a {T) ~ 7.7 x 10- 2 m a (T = 0)(A QC D/T) 3 - 7 
Axion oscillation starts at T a for which m a (T a ) = 3H(T a ): 

•" U)0 MeV J UoOGevJ F„ 



We refer the reader to paper |18| for the available formulae. If T a > Trh-, the coherent axion 



energy density is diluted by the entropy produced between T a and T RH . At the end of the 
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entropy dumping (around Trh), the coherent axion number density in unit of the entropy 
density is given by Yf ~ 6 2 F a 2 m a (T a )R a '/ Sf where 9 denotes the initial misalignment angle of 
the axion field, R a is the scale factor at T a and Sf is the total entropy at Trh- The ratio of 
the axion energy density to the critical energy density at present is given by 

/ w \ L5 / r \ °- 98 /A \ - 1 - 9 

~ 1 f lO^eV y- 44 f^^ 9 f Ago, y 19 (ig) 

V^J V F « / V300GeVy V200 MeW 1 ; 

where we have used ~ B~ 1 M^/32irF a 2 . The above result is valid only for n > 2 yielding 
Trh <T a . As we have anticipated, it shows that the case of n = 2 or 3 with F a ~ (M ¥ ,Mp) 1 / n+1 
yields a coherent axion energy density not exceeding the critical density although the corre- 
sponding F a exceeds 10 12 GeV. Furthermore, in this case of n = 2 or 3, axions can be a good 
dark matter candidate for an appropriate value of M v , which was not possible for n — 1. 

We remark that diluting the coherent axions with Trh < T a is allowed only when R-parity 
is broken. If not, stable lightest supersymmetric particles (LSP) produced after the flaton 
decay would overdose the universe. This can be avoided if the reheat temperature is bigger 
than the decoupling temperature of LSP which is typically Mlsp/20. However this is usually 
above 1 GeV, i.e. above T a . Consequently, the usual upper bound F a ~ 10 12 GeV can not be 
relaxed when the reheat temperature is bigger than Mlsp/20. We stress here that even when 
R-parity is broken and thus LSP cannot be a dark matter candidate, coherent axions can be 
a viable dark matter candidate when n = 2 or 3. 

Baryogenesis. Thermal inflation driven by PQ flatons may dilute away any pre-existing 
baryon asymmetry. However, PQ flatons themselves can produce baryon asymmetry after 
the reheating through the DH mechanism ||. A complicated Affleck-Dine type baryogenesis 
after thermal inflation has also been explored in Ref. ||19|| . Our previous discussion of flaton 
couplings in DFSZ type models indicates that flatons going to stops can be the most efficient 
decay channel. The decay-produced stops subsequently decay to generate a baryon asymmetry 
provided that the baryon- number violating operator, e.g., X'^U^D^D^ and the corresponding 
complex trilinear soft-term are present. Note that the PQ symmetry [see Eq. (0)] can be 
arranged so that dangerous lepton-number violating operators LQD C , LLE C are forbidden for 
the proton stability. 
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In order for the baryon asymmetry not to be erased the reheat temperature ([5]) has again 
to be less than few GeV 0. This again means that the DH mechanism can work only for 
n = 2 or 3 [see Eqs. (3) and (5)]. The produced baryon asymmetry is 

^^ 3 7W A (2Q) 

where AB is the baryon asymmetry generated by each flaton decay into stop-antistop pair. 
Using Eq. ([5]) and the estimate of AB given in |§, we find 

l A 332l 1n _ 2 d 17 Unn^y ' 



3xlo -io i 33.1 ^ 10 _ 2 ) \B a J \ F a J V 300 GeV, 
where arg(Am^ 2 ) denotes the CP violating relative phase which is constrained to be less than 



10 for superparticle masses of order 100 GeV [20|. For n = 3, the desired amount of baryon 



asymmetry can be achieved only when A332 is of order unity, while for n — 2 it can be done 
with a smaller X'^ 32 - 

In conclusion, we have examined some cosmological consequences of supersymmetric axion 
models in which the axion scale is radiatively generated as F a ~ (ma^MjS) 1 /™" 1 " 1 . In such mod- 
els, the early universe inevitably experiences a period dominated by the coherent oscillation of 
PQ flatons which start to oscillate at temperature around m 3 / 2 . Then a significant amount of 
oscillating PQ flatons can decay into axions, thereby yielding a too large axion energy density 
at the time nucleosynthesis. This consideration puts a limit on the effective branching ratio 
B a measuring how large fraction of oscillating flatons are converted into axions: it should 
be less than 1/3 ~ 0.02 depending upon our choice of the allowed extra number of neutrino 
species 5N U = 0.1 ~ 1.5. Models of hadronic axion with a radiative mechanism would yield 
B a ~ 0.5 unless the flaton coupling to the lightest flatino is unusually large. This is essentially 
because the flaton couplings to SSM are loop suppressed compared to the couplings to PQ 
sector. DFSZ type models with a radiative mechanism is more interesting since it can provide 
a rationale for the size of the /i term (and also the scale for neutrino masses). If the flaton 
mass M 9 2> M SSM denoting the typical mass in supersymmetric standard model, DFSZ type 
models would also suffer from the same difficulty as that of hadronic axion models. However 
for My comparable to M SSM , requiring B a to be about 1/10 does not provide any meaningful 
constraint on DFSZ type models. If one wishes to achieve a smaller B a , say about 10~ 2 in 
DFSZ type models, one then needs a kind of tuning of the model. Flaton decays into the 
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lighter stops is then picked out as one of the efficient decay channels leading to such a small 
B a provided 4M? < M 2 < \M\. 

Another interesting cosmological consequence of decaying flatons is the relaxation of the 
cosmological upper bound on the axion scale. For the axion scale bigger than 10 12 GeV, 
the entropy production by PQ flatons ends after the axion field starts to oscillate by QCD 
instanton effects, thereby diluting the coherent axion energy density in a rather natural way. 
With this late time entropy production by PQ flatons, the upper bound on the axion scale F a 
can be pushed up to about 10 15 GeV, but at the expense of breaking R-parity to avoid a too 
large mass density of relic LSP. Then the integer n which determines the axion scale in terms 
of m 3 / 2 and M P can take n — 1, 2 or 3. 

It is likely that any pre-existing baryon asymmetry is completely diluted by the huge 
entropy dumping in thermal inflation scenario. As the PQ flatons are expected to decay 
dominantly into stops, the DH mechanism for the late time baryogenesis can work in a natural 
manner when n = 2 or 3 so that the reheat temperature does not exceed 1 GeV. With 
broken R-parity, LSP is no more stable and can not be a dark matter candidate. In this 
scenario, coherent axions can provide a critical mass density of the universe by saturating the 
cosmological bound on F a which now can be as large as 10 15 GeV. 

Interestingly enough, we now observe that the case of n = 2 or n = 3 provides a very 
concordant cosmological scenario: (i) a proper baryon asymmetry is generated by the DH 
mechanism using baryon-number violating interaction \"U C D C D C , (ii) potentially dangerous 
coherent axions (with F a ^> 10 12 GeV) are diluted by the late time entropy production, (iii) 
both the baryogenesis and axion dilution require R-parity to be broken, and then diluted 
coherent axions constitute dark matter in the universe, 
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